Double-Slit Interferometry with a Bose-Einstein Condensate 
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A Bose-Einstein "double-slit" interferometer has been recently realized experimentally by Y. 
Shin et. al., Phys. Rev. Lett. 92 50405 (2004)- We analyze the interferometric steps by solving 
numerically the time-dependent Gross-Pitaevski equation in three-dimensional space. We focus on 
the adiabaticity time scales of the problem and on the creation of spurious collective excitations as a 
possible source of the strong dephasing observed experimentally. The role of quantum fluctuations 
is discussed. 

PACS numbers: 



Introduction. Several current efforts in the field of 
dilute Bose-Einstein condensates (BEC) are focusing on 
the creation of new technological devices, including quan- 
tum computers (jj and ultrasensitive interferometers Q 
to detect and measure weak forces. Atom wave interfer- 
ometry already provides unprecedented sensitivities to 
detect rotations, accelerations, and gravity gradients 
Performances can be further improved with interferome- 
ters based on BEC, which are the highest brilliant coher- 
ent sources of matter waves and which allow for larger 
separations between different interferometric paths. 

Interference between two spatially separated conden- 
sates was first demonstrated in , and theoretically an- 
alyzed in 7]. A BEC trapped in a harmonic magnetic 
trap was split in two symmetric halves by a laser knife. 
After releasing the external fields, the interference pat- 
tern of the overlapping condensates was observed with 
destructive imaging. The measured relative phase, how- 
ever, was not reproducible from shot to shot but was ran- 
domly distributed due to the presence of a large noise in 
the relative positions of the laser and the magnetic trap 
and of parasite currents when switching off the magnetic 
fields. Reproducible interference patterns were eventu- 
ally observed by trapping a condensate in deep optical 
periodic potentials ^ . Neighboring wells of optical lat- 
tices, however, are separated only by a fraction of a mi- 
cron and cannot be individually addressed, limiting their 
applications in technological devices. 

Recently, a stable double-well trap has been created 
in 0. A single coUimated laser beam was split with 
an acoustic-optical modulator, and finally focused by a 
lens. A single, cigar-shaped condensate was trapped in a 
double-well having a barrier much smaller than the BEC 
chemical potential, see Fig. (^). The condensate was 
then split along the axial direction by linearly increas- 
ing, in a ramping time tramp, both the distance between 
the two wells and the height of the interwell barrier, see 
Fig. CJd). The final distance between the two conden- 
sates was 12 /im, allowing for individual addressing 
and manipulation. After holding the two condensates in 
the respective traps for a time thoid, the confining field 
was turned off. The interference pattern of the two over- 
lapping condensates was measured by destructive tech- 



niques and was reproducible in different realizations of 
the experiment. However, a loss of coherence was ob- 
served when the condensates were held in the separated 
wells longer than tuoid ~ 5 ms. This strong dephasing 
has been tentatively attributed to axial and breath- 
ing mode excitations created during the splitting of the 
condensate. Attempts to increase the adiabaticity of the 
process with larger separation times tramp did not im- 
prove the stability of the measured phase. The dephas- 
ing manifested itself as a gradual decrease of the phase 
contrast accompanied by a bending and kinking of the 
interference pattern, which eventually makes the phase 
measurement impossible. 

In this Letter we theoretically analyze the MIT ex- 
periment as a prototype of a general BEC interferome- 
ter. Our analisys is quite general, and is relevant, for 
instance, in the study of the splitting and recombina- 
tion of BEC propagating in atom-chip wave guides j^]- 
We focus on the role of the collective excitations created 
during the splitting process as a possible dephasing mech- 
anism. We numerically solve the full three-dimensional 
time-dependent Gross-Pitaevskii (CP) equation to repro- 
duce the interferometric steps realized in the experiment. 
We study the adiabaticity of the splitting process and 
predict the interference contrast as a function of the var- 
ious time scales of the problem. We finally consider the 
corrections to the GP dynamics arising from quantum 
fluctuations. 

Protocol. Our simulation of the interferometric ex- 
periment consists of four steps: i) Initialization. We load 
a BEC into an optical double well potential of gaussian 
shape in the x direction and of harmonic shape in the 
y and z directions with frequencies ujy and respec- 
tively. We fix an initial separation xq between the wells 
in such a way that the height of the potential barrier is 
much smaller than the chemical potential, Fig. (^). We 
solve the GP equation in imaginary time in order to find 
the ground state (GS) of the system, ii) Separation. We 
separate the wells by ramping linearly from at time 
t = to Xramp at tramp, Fig. 1^). hi) Holding. Once 
the wells are separated, we generally allow the wavefunc- 
tion to evolve in the trap for a time thoid- iv) Ballistic 
(free) propagation. After a time t = tramp + thoid, we sim- 
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FIG. 1: Double Well potential in the x direction, a) The ini- 
tial condensate with a distance between the two wells of xq = 3 
fim and the chemical potential higher than the barrier height, 
b) Displaced wells in a time t — tramp with xo{tram,p) = 6 ^m. 
In this configuration we have two independent condensates. 



ply release the trapping potential and allow the packets 
from the separated wells to merge and overlap, generat- 
ing an interference pattern. The condensates are finally 
imaged after a time tfree of ballistic expansion by inte- 
grating along the y direction, simulating in this way the 
data collected experimentally. We cast the interference 
problem into a three-dimensional Gross-Pitaevskii (GP) 
equation. We introduce scaled units with the energy in 

-, and time in ujz^^, obtaining 



fkOz, length in dz 



h 

muJz 



-iv2 + K.t(r,t) + VVL(r,i) (1) 



where ViVL('",i) = 47rA^j-|?/'(T",t)|^ is the nonhnear po- 
tential arising from the interatomic interaction, with a 
the scattering length and N the number of condensate 
atoms in the trap. The external potential Vext{f,t) 
is given by a combination of a harmonic confinement 



along the y, z-directions, Vho{y, 



and a double-well time-dependent gaussian confinement 
along the x direction Vd^{x, t) = --^ ^-ix-xo{t)f /2a^ _^ 

g-(a:+a;o(t))^/2(T^ ^ whcrc 2xo{t) is the distance of the po- 
tential wells at time t, and cr is the width of each gaussian 
well. We base our simulations on the MIT experiment 
considering TV ^ 510^ ^'^Na atoms. The interatomic scat- 
tering length is a = 2.8 nm, and the trap parameters are 
ujy = 615 X 27r Hz, (jj^ = 30 x 27r Hz (giving a length unit 
dz — 3.8 /im), Vq ^ h X 5 kHz, and a = 2.5 /im. The 
bottom portion of each gaussian well approximates a har- 
monic potential with frequency cj^ = ^ ma'^' ~ ^ 

593 

Hz In the initial configuration, Fig.(^a), xq = 3 [im and 
the chemical potential fi = h x 1.78 kHz. The wells are 
then separated to a distance xo — 6 /im. Fig. (Qb), 
creating two condensates each with chemical potential 
H = hx 2.04 kHz. 

Results and Discussion. We performed GP simula- 
tions by varying the parameters tramp and t^nid keeping 



fixed the expansion time i/r 



5.3 ms 



tral goal of this section is to extract the relevant adia- 
baticity time scale of the double-slit BEG interferometer. 
This requires studying the excitations created during the 
splitting of the trap and how these propagate during the 
holding time. 

Ramp time. We first examine the behavior of the sys- 
tem as a function of the separation time, t 
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for fixed 

thoid = 10.6 ms. In Fig. ((2Jl, we show the contours of the 
y-integrated probability density Py{x, z, t), which closely 
resemble the projected images of the experiment. For 
short separation times we have observed considerable dis- 
tortions and dephasing in the interference patterns. 

Given the symmetries of the system, low-lying excita- 
tions can only have even parity. For example, in the axial 
direction the longest oscillation period is ^ = 16.6 
ms, and faster ramping times can easily excite large am- 
plitude monopole and quadrupole oscillations. However, 
we have numerically checked that such collective modes 
retain a high degree of collectivity (due to the harmonic- 
ity of the trap) and cannot be a source of dephasing. 
Nonlinear couplings, which can potentially destroy the 
interference patterns, apparently occur at much longer 
times than those considered in the experiment. 

We now consider the frequencies of low-lying even ex- 
citations along the radial direction as a function of the 
interwell distances. The oscillation periods first increase 
along with the increase of the interwell barrier height. 
These reach a maximum value when the difference be- 
tween the energy of the interwell barrier and the chem- 
ical potential becomes equal to the frequency of the 
first odd collective mode, corresponding to the Joseph- 
son "plasma" oscillation. Further increasing the interwell 
distance results in the periods decreasing to a final satu- 
ration once the two wells are completely separated. This 
behaviour reflects the evolution of the single particle en- 
ergies of the double well during the ramping process. The 
maximum oscillation period sets the relevant adiabatic- 
ity time scale of the double-slit BEG interferometer. 

In each separate well, both even and odd excitations 
modes can be excited along the radial axis. The low- 
est in energy is the dipole collective mode, whose oscil- 
lation frequency can be calculated in a variational ap- 
proach to be ujjy = Lo^ { i+^\/cr2 )^^^'^\ with 7 the the 
radial width of the condensate. The trapping poten- 
tial is highly anharmonic and the correponding oscilla- 
tions frequencies strongly depend on the chemical po- 
tential of the system, even in the Thomas-Fermi limit 
(contrary to what happens with an harmonic confine- 
ment). With N — 10''', the largest frequency is obtained 
at an interwell distances xq ^ 3.7 /im and chemical po- 
tential /i ~ 0.64 Vdw{x = 0): uimin = 0.65 w^, and 
the oscillation period is r^ax — 1-54 r^. For compari- 
son, when Xq = 6 ^m, and the two wells are completely 
separated, ujd = 0.75 and the oscillation period is 
TD = 1.33 = 2.58ms. 

The breaking down of adiabaticity, however, is not the 
only source of dephasing. We have not observed any loss 
of visibility after imposing harmonic confinement along 
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FIG. 2: Contour plots of the probability density integrated 
over the y-coordinate Py{x,z,t) from 3D simulations for 
thoid= 10.6 ms and t/ree ~ 5.3 ms, as a function of the ramp- 
ing time tramp'- a) 2.7 ms, b) 5.3 ms, c) 10.6 ms, and d) 21.2 
ms. 



FIG. 3: Contour plots of the probability density integrated 
over the y-coordinate Py{x,z,t) for tramp= 5.3 ms and t/ree 
= 5.3 ms as a function of the hold time thoid- a) 0, b) 5.3 ms, 
c) 10.6 ms, and d) 21.2 ms. 



the X axis. On the other hand, we have found that, in 
the presence of the gaussian confinement, the interference 
pattern quickly degrades while increasing the number of 
atoms, keeping all other parameters fixed. Fast ramp- 
ings induce large amplitude oscillations along the radial 
direction of each trap. Due to the combined anharmonic- 
ity of the confining potential and the non-linearity arising 
from the interatomic interaction, the energy of the excita- 
tion quickly redistributes among different modes, rapidly 
damping the collective oscillation. We have not observed 
any coupling of such modes with oscillations among the 
radial axis. It is useful to recall that the total energy 
of the system remains conserved in the GP simulation, 
and the effective damping is only due to redistribution of 
the collective energy among different Bogoliubov excita- 
tions levels. With smaller oscillation amplitudes, induced 
with longer ramping times, the visibility of the fringes 
has been greatly improved Fig.®. However, this effect 
is in contrast with the experimental results where an 
increasing of the ramping time, even well beyond Tmax, 
does not improve the interference pattern, which degra- 
dates for thoid > 5 ms (see Fig. 3 in 01), independently 

Hold Time. We now investigate the role of the hold 
times thoid in the deterioration of the interference pat- 
tern. Fig. (PI presents contours for the y-integrated prob- 
ability density for different thoid and fixed tramp — 5.3 
ms. We clearly note a growing deformity of the patterns 
as the confinement time increases, along with the redis- 
tribution of the initial energy among various modes. The 
maximum dephasing indeed occurs when the oscillation is 
completely damped. Our simulation goes in the direction 



of the experiment: an increasing of t^oid corresponds to a 
growing distortion of the fringes pattern. However, there 
is an important difference. We do not observe a complete 
distortion of the fringes pattern for thoid > 5 ms, as in 
01 . Instead, we find distorted but still defined fringes 
even for thoid ~ 20 ms. Moreover, such distortion can 
eventually be reduced by increasing the ramping time, 
as suggested by Fig.©. Changing the relative depth of 
the two wells to include an asymmetry in the x-direction, 
has not lead to a more rapid deterioration of the fringe 
contrast. 

Another important difference between the experimen- 
tal findings and our GP simulations occurs in the shape of 
the distorted interference pattern. We always find sym- 
metric bending or kinking of the fringes, contrary to that 
observed in Fig. (3) of jjj. The asymmetric kinking of 
the experiment suggest the presence of slight additional 
geometry variations of the two traps. In the real exper- 
iment, indeed, the shape as well as the parameters of 
the trap can be slightly different from the corresponding 
mathematical form we have used in our simulations. A 
small disalignement of the two traps, for instance, can 
induce small amplitude out of phase dipole oscilations 
along the z axis which can originate an asymmetric bend- 
ing. Residual excitations due to the transportation of the 
condensate to the "science chamber" Q prior the initial- 
ization step i) , might also be the source of the dephasing 
observed experimentally, (l^j ■ We can conclude this part 
by stating that the excitations coming from the creation 
of two separated condensates give significant distortion 
and loss of contrast in the interference fringes, but they 
do not destroy completely the pattern. In general such 
distortions increase with thokh as in the experiment, and 
with the chemical potential, but it can be substantially 
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reduced by decreasing the ramping time tramp- 

Beyond Gross-Pitaevskii. The Gross-Pitaevskii 
framework does not account for dcphasing arising from 
quantum (many-body) correlations. While, on the one 
hand, a full calculation of the quantum dynamics is well 
beyond current computational capabilities, simple esti- 
mates [iSl [l^ [iSl [IGJ of the dynamical evolution of the 
fringe contrast can be carried out within the two-mode 
approximation and tested experimentally. The quantum 
phase Hamiltonian governing the condensate is (setting 
n= 1) [13: 
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where 'i>{(j),t) is the probability amplitude for the rela- 
tive phase (p of the two condensates. The "charging" en- 
ergy Ec = 2^ remains approximately constant during 
the splitting of the condensate, while the Josephson cou- 
pling energy Ej decreases exponentially with the inter- 
well distance. In the strong coupling limit Ej/Ec >> 1 
(achieved when the chemical potential is close enough 
to the interwell barrier), the phase oscillation frequency 

4£'|/iV2. The phase probabihty 

|4'((/), t)!'^ has a gaussian distribution with dispersion 
— << I. As long as a{t) remains small dur- 

ing the dynamics, the system can be described in the GP 
framework. With a linear ramping, d{t) = 2xot/Tramp, 
and in the WKB approximation, the Josephson coupling 
energy is Ej{t) ^ e("*/'re//)^ with r^ff = Tramp/ S and ef- 
fective action S = 2a;o\/2m(Vb — /j,). At a first stage, 
while Ej decreses with time, the amplitude 5'(^, t) fol- 
lows adiabatically the ground state of the effective quan- 
tum potential in Eq.|2Jl. Breakdown of adiabaticity oc- 
curs at the freeze-out time tj given by ujj{tf) ~ TJJ^ 
namely when the Josephson period becomes equal to 
the characteristic time of the change of the Josephson 
coupling energy Te/f = — (^^/i?j)~^. At longer times 
the dynamics is dominated by the kinetic part of the 
quantum Hamiltonian Eq.lj^J. The temporal evolution 
of the phase fluctuations is therefore given by cr^(t) = 
a'^{tf) + El/'i<T'^{tf)t^, with (j{tf) being the width calcu- 
lated at the freeze-out time. Replacing the corresponding 
expressions, we finally obtain: 



a'{t) = |^(c.,.(i/)-^+-.(i/)^^) - f 



^2 



' ramp 

(3) 

which is the central result of this section. To experimen- 
tally test the quantum phase dynamics Eq.lPJ it would 
be necessary to average over several identical interfero- 
metric realizations. In each experiment, an interference 
pattern will actually be observed 0], but with a rela- 
tive phase chosen randomly with a gaussian distribution 
of width given by Eq.Q. The quantum dynamics will 
therefore manifest itself as a loss of fringe visibility of 



the ensemble averaged interference patterns, which will 
be almost completely washed out when a ^ 2. Replac- 
ing the experimental values of 4], this happens with a 

holding time thoid - \J ■:^§^Jramp — 27r50 ms. 

Conclusions. We have theoretically studied the 
double-slit interferometer recently created in ,4j solving 
numerically the full 3D time-dependent Gross-Pitaevski 
equation. We have studied the adiabaticity time scales of 
the system, concluding that the dephasing arising from 
the creation of spurious excitations can be strongly re- 
duced by increasing the ramping time of the double- well. 
Such findings indicate that the loss of visibility and the 
bending of the interference pattern observed experimen- 
tally should arise from a different noise source. 
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